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Abstract 

Maximal Abelian gauge fixing and subsequent Abelian projection 
of SU(2) lattice gauge theory defines closed trajectories of magnetic 
monopoles. These trajectories can be interpreted in terms of an effec- 
tive scalar field theory of the MAG monopoles using the worldline rep- 
resentation of the functional determinants. Employing the monopole 
worldlines detected in the numerical simulation, we show that a scalar 
bound state exists. The screening mass m of this state properly scales 
towards the continuum limit. We find m ~ 1.3 GeV when the string 
tension yfa = 440 MeV is used as reference scale. 
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Introduction. 

Lattice calculations |]-@ performed in the so-called "Abelian gauges" || |J 
have provided evidence that a condensate of magnetic monopoles exist in the 
Yang-Mills vacuum. Consequently, the vacuum expels color-electric flux by 
virtues of the dual Meissner effect and produces confinement. In this sense, 
the vacuum represents a dual superconductor. On a phenomenological level 
a superconductor can be described by a Ginzburg Landau theory. There 
have been attempts to construct the pertinent dual Ginzburg Landau theory 
for the QCD vacuum or to extract it from lattice gauge simulations [|T0||— [|T3[ . 



The difficulty seems to consist in the mapping of the monopole degrees of 
freedom to the one of the scalar Ginzburg Landau field. 

The Ginzburg Landau theory describes a complex scalar field interacting with 
an electromagnetic gauge field. The (dual) Abelian electromagnetic field can 
be, in principle, be integrated out yielding an effective theory of a complex 
self-interacting scalar field. In this paper, we will extract the effective scalar 
field theory describing the dual superconductor of the QCD vacuum with 
the help of lattice gauge simulations. To this end we firstly determine the 
ensemble of magnetic monopole loops of the vacuum by performing a lattice 
calculation in the maximum Abelian gauge, performing the Abelian projec- 
tion and identifying the magnetic monopole loops by the method of DeGrand 
and Toussaint ]14|]. The obtained ensemble of closed (magnetic monopole) 



trajectories is then described in terms of an effective scalar field theory by 



using the worldline formalism [15, 16]. Our lattice simulations will show that 



a scalar anti-scalar bound state survives in the continuum limit. 

The Ginzburg Landau theory of MAG monopole trajectories. 

The central idea of the present paper is that the theory of closed monopole 
trajectories arising from MAG projected SU(2) lattice gauge theory is equiv- 
alent to a theory of a charged scalar field. This scalar field theory designed 
to describe the monopole properties of the Abelian projected SU(2) Yang- 
Mills theory necessarily inherits the scaling laws from the underlying SU(2) 
Yang-Mills theory and, in particular, the property of asymptotic freedom by 
construction (provided that SU(2) monopole theory properly scales towards 



the continuum limit). It was pointed out by Zakharov [|18|] that the scalar 
field theory which emerges from SU(2) monopole loops is an interesting can- 
didate for avoiding the so-called fine tuning problem, which is generic in 
(4-dimensional) scalar field theory equipped with the standard 4> 4 potential. 
Since scalar field theories with local interactions of the scalar field possess 
an infra-red fixed point, the action term is presumably not a polynomial of 
finite order. One might argue that the increase of complexity due to the 
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non-local interactions prohibit the access to such theories at a practical level 
(e.g., the numerical simulation). However, examples of scalar theories incor- 
porating asymptotic freedom at the expense of the non-locality of the action 



have been treated in the literature 19 - 21 



In order to establish the equivalence between the theory of the monopole 
loops and the scalar field theory, we consider the general form of the partition 
function of a complex scalar field 



Z[M] = J V<j>V<fiexp{ - I d 4 x$(x) -d 2 + m 2 + M(x) <j>{x) (1) 



Here m is the usual mass term, and M(x) is an external source, which we 
will specify later. V($4>) describes the interaction of the scalar field. The 
only restriction which we impose here is that we assume the potential term 
V(<j)'<j)) to admit a Taylor expansion so that ([!]) can be written as 



Z[M) 



cxp 



d A x V( 



S 



8M(x) 



Z [M] , 



(2) 
(3) 



Z [M] = Det _1 [-9 2 + m 2 + M(x)} 
Using the proper-time representation of the functional determinant in @ 



i.e., 



r [M] 



\nZ [M] = J ^e~ m2T tr exp{-r(-<9 2 + M(x) ) } , (4) 



the emerging heat kernel can be interpreted as the time evolution operator 
of a point particle, for which the usual Feynman path integral representation 
holds (we refer to |E2| for a recent review of the world line formalism) 



r [M] 



dT 



-m 2 T 



/ d 4 x / Vx{t) e o v ■ (5) 



x(T)=x(0) 



Here we have split off the integral over the zero-modes of the path integral, 
/ d A xo, where xo, the so-called loop center of mass, corresponds to the average 
position of the loop: Xq := (1/T) Jq dr x M (r), i.e., 



Vx{t) 



Hdx^r) 5 (4) [x^ - (1/T) / drx^r) 
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In order to relate the functional integral over the world lines x(r) in (|5|) to 
the expectation values over loop clouds, we normalize it with respect to the 
free theory (M = 0) and introduce 

t .2 
— f dr — 

(0(x)) x = AT 1 J Vx(r)e o 4 o(x{t)) . (6) 

x(T)=x(0) 

where 



T ,2 



. * t . , ~I dT2 t r d 4 p n 2 T 1 

Equation (|^) defines the expectation value of an observable O evaluated over 
an ensemble of closed loops x(r); the loops are centered at a common average 
position xo ( "center of mass" ) and are distributed according to the Gaussian 
weight exp[— f dr^A. These definitions lead us to the compact formula 

oo 

r [M] = t^t 2 J d'x J ^ e- m2T (exp{- J M(x(r)) rfr} 



o 



The world line representation of the interacting scalar field theory (0) is 
obtained by inserting (^) into (0). Thereby, the interaction of the scalar field 
V(4>^4>) gives rise to an effective interaction of the loops V(x(t)), i.e., 



oo 

T[M] = j±-J d*x |-e" m2T (e-v<*r» ex p{"/ M ^)) & 

(9) 

Note, that once a particular choice of the scalar interactions V(0V) is made, 
the determination of V(x(t)) in (^) is in principle straightforward. 

In this paper, we propose to determine the closed monopole loops of the 
Yang-Mills vacuum using lattice gauge simulations (see below for details). 
Employing the above illustrated equivalence between a world line ensemble 
and a scalar field theory, the effective scalar field theory underlying the dual 
superconductor of the Yang-Mills vacuum can be in principle extracted as 
suggested in [|18|. It was recently observed that Monte-Carlo calculations of 
the loop averages (such as those in @ and (||)) are feasible [|23], |24 . 



In order to connect properties of the monopole loops with expectation values 
of the scalar field theory, we study different choices of the external current 
M(x) in the context of the scalar field theory (0H) and in the context of the 
worldline formalism (P), respectively. 
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To shorten the presentation, we introduce the shorthand notation 

OO 

0( 4 = wfl A » / % *~ m * T ( e " WT)) ( i(t) )), • (10) 

Firstly, we choose 

M(x) = j5\x-x ) (11) 
and insert this ansatz into ([D yielding 

dr M d *zw = (M*,)) ■ (12) 



dj dj 

On the other hand, it is clear from (|9]) that dT[M]/dj counts the number 
times a monopole loop passes through the specified point Xo, 

p(x ) = (^J dr6 4 (x(r)-x ) , (13) 

and, hence, corresponds to the probability of finding a monopole or anti- 
monopole (depending on the orientation of the trajectory) at xq. We will 



call this quantity monopole density. Comparing (|T^) and (|12|), one identifies 

p(x ) = ^(x )cf>{x )^ . (14) 

In order to get a first insight into the propagators of the full interacting scalar 
theory (|l|), we investigate the particular choice of the source 

M(x) = ji 5\x - x ) + 32 S\x - y ) . (15) 

Inserting (|15D into (|l]), taking the derivative with respect to the currents j\ 
and 32, respectively, we obtain the connected Green function 

C(x -y ) = -f—\nZ[M] (16) 
dji d3 2 

^<f>{x ) 0V(l/o)\ - (tf<t>(xo)) U j <f>{yo) 



which in view of ( |T3"D can be interpreted as the correlation function of the 
monopole density. The long distance behavior of this correlation function 
determines the screening mass of the scalar anti-scalar excitation. 
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Figure 1: Closed world line contributing to the scalar correlation function 
C(x — y ) (left panel). The quark anti-quark potential as function of the 
quark anti-quark distance r (right panel). 



Taking the derivatives with respect to the currents ji and j'2 of the effective 
action @ in the world line formulation yields the loop representation of the 
above correlation function: 

C(xo-yo) = (^Jdr5\x(T)-x Q ) J dr' 5 4 {x{r')-y Q )Jj (17) 

The latter equation has a simple interpretation: the correlation function 
C(xq — yo) is obtained by taking the average over all closed loops which pass 
through points xq and yo, respectively (see figure [l|). Once the monopole 
world lines are at our disposal, we are able to calculate the full propagator 
( ]T6l) of the corresponding scalar field theory without specifying the scalar 
interactions V{(jr 0). 

Lattice results. 

In order to determine the closed worldlines of the magnetic monopoles we 
performed simulations on a 12 3 x 24 lattice using the Wilson action. In order 
to express the size of the lattice spacing in physical units, we will make use 
of the 1-loop scaling relation 

aa 2 ((3) = s expj-^ (/? - 2.3) J , (18) 

where So is the parameter expressing the string tension in units of the lattice 
spacing for (3 = 2.3. The scaling relation is assumed to reproduce the string 
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tension within the scaling window (3 G [2.1,2.6]. Since we are using an 
asymmetric lattice, we have checked this assumption by calculating the quark 
anti-quark potential V(r) for the above lattice size. For this purpose, we 
calculated the Polyakov loop correlation function P(r) as function of the 
quark ant i- quark distance r, i.e., 

P(r) oc exp{- 24 V(r) a(P)\ . (19) 



For this task, we used the 2-level Luscher-Weisz method [2S|. The averages 
at level 1 were performed using 50 iterations while 10 evaluations were em- 
ployed for the averages at level 2. 150 independent 2-level measurements 
were performed to yield the accuracy of the data shown in figure p] (right 
panel). The raw data were fitted to the potential form 

V(r) = -- + ar . (20) 
r 



Fit parameters have been the parameter so in (|18| ) and the offset of the 
potential for each {3 value. We finally obtained: 

s = 0.146 ±0.005 , a = 0.183 ±0.005 . (21) 

This value for s is in reasonable agreement with the known value 0.125 for 
symmetric lattices The agreement between the lattice data and the 



model potential (|20"D is very good (see figure |l|). In the following, we will use 
a string tension a = 440 MeV as reference scale. 

The Maximal Abelian gauge (MAG) condition, i.e. 

J2 trfc(x) r 3 U n ^(x) r 3 ) $ max , (22) 
{x},ti 1 ' J 



where U^(x) = Q(x) U^x) fl^x+fi) are the gauge transformed link variables, 
is implemented by employing a standard iteration over-relaxation algorithm. 
We do not expect that our numerical procedure locates the global maximum 
of the non-linear functional (|22|) . Choosing different sets of local maxima of 



(|22|) implies that different gauge conditions are implemented (see e.g. |27| for 
a more detailed study of the issue gauge fixing ambiguities). We stress that 
the properties of the monopoles corresponding to these different gauges might 
turn out to be different. Rather than pursuing a detailed study of the effects 
of these so-called Gribov ambiguities, the aim of the present paper is to show 
that a monopole anti-monopole bound state exists at least for a specific choice 
of the gauge (MAG, using an iteration over- relaxation algorithm to install the 
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gauge condition). Once the MAG is implemented, the SU(2) gauge theory is 
projected onto a compact U(l) gauge theory by the usual prescription 

U%(x) = exp{^V a } -> exp{^ 3 r 3 } . (23) 

Once the compact U(l) gauge theory is at our disposal, we use the stan- 
dard method of DeGrand and Toussaint [HJ to extract the closed monopole 
trajectories. 

In order to obtain the screening mass m of the scalar bound state, instead 
of ( |I7| ) we consider the related correlation function 

C(t) = (f(t) m) - (m) (m) « exp{-mi} , (24) 

where 

m := "£tf(t,x)<f>(t,x). (25) 

{£} 

The correlation function C(t) is obtained from the monopole trajectories as 
follows: the number n t of monopoles is counted for a given time slice. The 
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correlation function 

C dis (t) = (ji t n \ (26) 

provides the disconnected counterpart of the Green function (|24j ) . The func- 
tion can be well represented by the fit function 

C dis (t) = p 2 + a exp{-mt} , (27) 

where we used the fact that the disconnected correlation function asymptot- 
ically (t ^> 1/m) approaches the monopole density squared. 

After thermalization, we performed 100 measurements which were separated 
by 15 dummy sweeps to reduce the auto-correlations. For each /3, we extract 
ma from a fit of ( |2"T| ) to the numerical data for C d i S {t) (Jffi). The figure ^ 
shows the normalized connected correlation function 

C n {t) = (c dis {t)-pA/a = expj-mt}, (28) 

where the value of t (in physical units), t = a{[3) N t , is obtained by using the 
scaling relation (|Tj|). We observe that the data points obtained from several 
(3 values fall on top of the same curve. This signals that the screening mass 
extrapolates to the continuum limit. Using a string tension o = (440 Me V) 2 , 
we find a mass m si 1.3 GeV, which is of order of the mass of the low lying 
glueballs. 

Acknowledgments: 

We thank V. I. Zakharov for stimulating discussions. 

References 

[1] A. S. Kronfeld, G. Schierholz and U. J. Wiese, Nucl. Phys. B 293, 
461 (1987). 

[2] T. Suzuki and I. Yotsuyanagi, Phys. Rev. D 42, 4257 (1990). 

[3] S. Hioki, S. Kitahara, S. Kiura, Y. Matsubara, O. Miyamura, S. Ohno 
and T. Suzuki, Phys. Lett. B 272, 326 (1991) [Erratum-ibid. B 281, 
416 (1992)]. 

[4] T. L. Ivanenko, A. V. Pochinsky and M. I. Polikarpov, Phys. Lett. 
B 252, 631 (1990). 



9 



[5] G. S. Bali, V. Bornyakov, M. Muller-Preussker and K. Schilling, 
Phys. Rev. D 54, 2863 (1996) ||arXiv:hep-lat/9603012|| . 



[6] A. Di Giacomo, B. Lucini, L. Montesi and G. Paffuti, 
Phys. Rev. D 61, 034503 (2000) ||arXiv:hep-lat/9906024 l 



Phys. Rev. D 61, 034504 (2000) [|arXiv:hep-lat/9906025 l 



[7] V. Bornyakov and M. Muller-Preussker, Nucl. Phys. Proc. Suppl. 
106, 646 (2002) ||arXiv:hep-lat/0110209|| . 

S. Mandelstam, Phys. Rept. 23, 245 (1976). 

[9] G. 't Hooft, Nucl. Phys. B 190, 455 (1981). 

[10] G. S. Bali, C. Schlichter and K. Schilling, Prog. Theor. Phys. Suppl. 
131, 645 (1998) ||arXiv:hep-lat/9802005 . 



[11] M. N. Chernodub, F. V. Gubarev, M. I. Polikarpov and V. I. Za- 
kharov, Nucl. Phys. B 600, 163 (2001) ||arXiv:hep-th/0010265| . 



[12] M. N. Chernodub, F. V. Gubarev, M. I. Polikarpov and V. I. Za- 
kharov, Nucl. Phys. Proc. Suppl. 102, 347 (2001). 

[13] Y. Koma, E. M. Ilgenfritz, H. Toki and T. Suzuki, Phys. Rev. D 64, 
011501 (2001) ||arXiv:hep-ph/0l03l62| . 

[14] T. A. DeGrand and D. Toussaint, Phys. Rev. D 22, 2478 (1980). 

[15] M. Stone and P. R. Thomas, Phys. Rev. Lett. 41, 351 (1978). 

[16] S. Samuel, Nucl. Phys. B 154, 62 (1979). 

[17] K. Halpern and K. Huang, Phys. Rev. Lett. 74, 3526 (1995) 
||arXiv:hep-th/9406l99|| . 



[18] V. I. Zakharov, |arXiv:hep-ph/ 0202040 



[19] K. Halpern and K. Huang, Phys. Rev. D 53, 3252 (1996) [|arXiv:hep- 
th/9510240 1. 



[20] K. Langfeld and H. Reinhardt, Mod. Phys. Lett. A 13, 2495 (1998) 
||arXiv:hep-ph/97022TT [. 



[21] R. F. Langbein, K. Langfeld, H. Reinhardt and L. von Smekal, Mod. 
Phys. Lett. A 11, 631 (1996) ||arXiv:hep-ph/9310335|| . 



10 



[22] C. Schubert, Phys. Rept. 355, 73 (2001) [ |arXiv:hep-th/0101036 |. 

[23] H. Gies and K. Langfeld, 

Nucl. Phys. B 613, 353 (2001) ||arXiv:hep-ph/0l02l85|| , 
Int. J. Mod. Phys. A 17, 966 (2002) ||arXiv:hcp-ph/0112198|| , 
K. Langfeld, L. Moyaerts and H. Gies, |arXiv:hep-th/0205304 



[24] M. G. Schmidt and I. O. Stamatescu, |arXiv:hep-lat/0201002 



[25] M. Luscher and P. Weisz, JHEP 0109, 010 (2001) ||arXiv:hep- 
lat/01080T4 . 



[26] C. Michael and M. Teper, Phys. Lett. B 199, 95 (1987). 

[27] L. Giusti, M. L. Paciello, C. Parrinello, S. Petrarca and B. Taglienti, 
Int. J. Mod. Phys. A 16, 3487 (2001) ||arXiv:hep-lat/0104012|| . 



K. Langfeld, H. Reinhardt and J. Gattnar, Nucl. Phys. B 621, 131 
(2002) ||arXiv:hep-ph/0107141|| . 



11 



SU(2), 12 3 x24 



m=1.5 GeV 




0.1 0.2 0.3 0.4 0.5 

t[fm] 



